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abstract 

The impulsive propellant reorientation process is modeled using the 
ECLIPSE Code. A brief description of the process and the computational model 
is presented. Code validation is documentated via comparison to experimen- 
tally derived data for small-scale tanks. Predictions of reorientation per- 
formance are presented for two tanks designed for use in flight experiments 
and for a proposed full-scale OTV tank. A new dimensionless parameter is 
developed to correlate reorientation performance in geometrically similar 
tanks. Its success is demonstrated. 


INTRODUCTION 

The ECLIPSE Code, (Energy Calculations for Liquid Propellants in a Space 
Environment), is being developed as one component of the reduced gravity fluid 
management technology program being conducted by the NASA Lewis Research Cen- 
ter (ref. 1). The long range goal of a general tool for computational model- 
ing of liquid propellant behavior in a reduced gravity environment is being 
pursued by stages with each stage corresponding to a problem of current inter- 
est to designers of advanced spacecraft. The ability of ECLIPSE to model jet 
induced mixing in propellant tanks (refs. 2 and 3) and ptopellant tank self- 
pressurization (ref. 4) has been documented. The focus of the work being 
reported in this paper is the modeling of liquid motion induced by a sudden 
change in the acceleration environment. 

During coast in Low-Earth-Orbit (LEO), liquid propellants can collect in 
various regions of the propellant tank due to atmospneric drag on the space- 
craft. The process of positioning the liquid over' the tank out'et by firing 
auxiliary thrusters is known as impulsive reorientation or settling. Since 
impulsive reorientation requires the expenditure of propellant, it >s impor- 
tant to optimize the process to minimize the associated propellant require- 
ments. If the thrust level is too low, the propellant may not reposition. If 
it is too high, a large geyser may form and vapor pockets may be trapped in 
the pool. Proper spacecraft design and operation requires a good understand- 
ing of the process and the parameteis which control it. 


•Work funded by NASA Grant NAG3-578. 


Small-scale experiments have been performed in the NASA Lewis Research Cen- 
ter drop tower and in the Zero-Gravity Facility to examine liquid motion 
induced by accelerations which model the reorientation process. The experi- 
ments were performed in transparent tanks and the fluid motion induced by 
these accelerations was recorded via high-speed photography. These experi- 
ments identified liquid-vapor interface shapes for a zero-g environment 
(ref. 5) and records of fluid motion induced in partially filled tanks due to 
an imposed acceleration were produced (ref. 6). Further studies conducted by 
I. Sumner (ref. 7) examined the energy expended to accomplish reorientation. 

A performance map for the reorientation process in these small-scale tanks was 
developed. This reference provides substantial detail for a number of acceler- 
ation levels and tank fillings and the results reported in this reference have 
served as the basis for code validation. 

Prediction of reorientation performance has relied on two tools. The sim- 
plest is a rigid body dynamics analysis which assumes the propellant pool to 
behave as a single solid body (ref. 8), Although this analysis is simple to 
perform, it is such a coarse assumption that large safety factors must be pro- 
vided when it is used for design purposes. The second tool is described by I. 

Sumner (ref. 7) and is based on an empirical analysis developed by J. Salzman 
(ref. 9). The computational procedure is an empirically based approach using 
a Weber number criteria to preclude geysering and results in calculation of a 
liquid leading edge velocity. I. Sumner (ref. 7) extended the analysis to 
include small geysers. Although predictions based on this method correlated 
well with reported experimental results, it does not start from first princi- 
ples and its application is limited to similar geometries and conditions. » 


NOMENCLATURE 

a acceleration 

Bo Bond number 

F vol ume-of-f 1 ui d function • 

Fr Frcude number I 

g body force per unit mass, typically gravitational ^ 

h length dimension, typically tank length j 

1 

h* nondimensional i zed length : 

p pressure 

R tank radius 

Re Reynolds number , 

Se settling number | 

t,T time 



X and y velocity components respectively 
velocity of liquid leading edge 

spatial coordinates, cartesian or axi symmetr i c (x in radial direction) 

fractional cell volume open to flow 

dynamic viscosity 

kinematic viscosity 

switching integer 

dens i ty 

surface tension 


ECLIPSE CODE 


The ECLIPSE Code is a descendant of the family of SOLA Code: 


:en at 


the Los Alamos National Laboratory. In particu’d'', the NA'^a-v/qfzD Code (ref. 
10) was used as the foundation upon which ECLIPSE is being ouilt. The base- 
line code solves the laminar hydrodynamic problem using the Volume-0f-F1uid, 
(VOF), algorithm to determine the location of the free surface. The computa- 
tional model uses a VOF function. F, to track the free surface, and a cell 
blockage function, 9, to model partial cell blockage. Equations expressing 
conservation of mass, the force-momentum balance, and the F-transport equation 
can be writ ten : 
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Although many features such as heat transfer and thermodynamic models have 
been added during the development of ECLIPSE, only minor modifications to the 
baseline code were required to study impulsive reor ientation . The routine 
which provides an initial low-g fluid configuration has been modified so that 
the liquid can be positioned at either end of the tank. A set of variables 
have been incorporated to impose a time dependent acceleration environment on 
the fluid in the tank. The ability to terminate execution based on a criteria 
for completion of reorientation has been incorporated into the code. Finally, 
additional output options have been added to enhance the tracking of variables 
of specific interest to the study of the reorientation process. 


CODE VALIDATION 

Six cases were selected from I. Sumner (ref. 7) to serve as verification 
that ECLIPSE accurately models the impulsive reorientation process. A summary 
of the test conditions is presented in table I. The fineness ratio, FR, is 
the ratio of tank axial length to tank radius. The percentage of the tank 
occupied by liquid is recorded in the column labeled FL. TCTFE istrichloro- 
tr i f 1 uoroethane . R is the tank radius, and the balance of the variables are 
self explanatory. These cases were selected to provide a range of geyser 
formation from large to nonobservable. 

A typical computational mesh used to model these cases is presented in fig- 
ure 1. The liquid is shown at the top of the tank in a zero-g configuration. 

The mesh has been refined near the maximum tank radius (i.e., tank wall) in 
order to assure a smooth transition of the flow from the barrel section into 
the head. Although all figures presented in this report show a fuM cross- 
section through the tank, the code solves the problem in cylindrical coordi- 
nates the settling acceleration and resulting flow is ax i symmetri c . Therefore, 
only half the number of cells depicted are required to perform the 
computations . 

A sequence of flow fields computed for the conditions specified for Test 5 
is displayed in figure 2 . A comparison between computational prediction and 
experimental observation of geyser tip location as a function of time is pre- 
sented in figure 3. ECLIPSE predicts foimation of a geyser with a maximum 

height of 3.2 cm whereas I. Sumner (ref. 7) reports a maximum geyser height of 
3,3 cm. ECLIPSE predicts dissipation of the geyser into the rising free sur- 
face at approximately 1.05 sec whereas I. Sumner (lef. 7) reports thisevent 
at approximately l.'O sec. A comparison between computational pr-edkticn and 
experimental observation of geyser tip location as a function of time for Test 

1 is shown in figure 4. A. Patag <ref. 12) in nis thesis on the work presented 

herein, presents comparisons between the data from '■he experiment listed in 
table I and the computational predictions. In gene'al, reasonable agreement 
was obtained with a single ex'eot'on. The parameter's specified for Test 3 com- 
bine to induce large leading edge veiscities as the iquid moves a‘Cng the 
tank walls. At geyser incept'on at tne bottom of tne tank, the free surface 
model in ECuIPSE fails to track soriestiy the f-ji'matkn of j pool and the ini- 
tial growth of the geyser. Comoii';ed hi)w '^ields for Test 3 are shewn in fig- 
ure 5. It should be noted that ’’he i-inetic erei'gy 'mpa-ted to the liquid is 
hign. As such, it 'S highly i ne*"’* ' c i ent and would tequi'e an excessive expen- 
diture of propellant to produce these conditions in a r ea ' spacecraft. Since 
the current study is focused on optimization of the - eor i entat ion process, the 
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inability of ECLIPSE to successfully niodel Test 8 is not viewed as a signifi- 
cant handicap. 

Based on the evidence presented in the preceding paragraphs, it was con- 
cluded that ECLIPSE is a suitable tool for modeling pulsed settling. 


SCALE MODEL OTV TANKS 

One component of the NASA Lewis reduced gravity fluid management technology 
program is development of a flight experiment to examine a variety of fluid 
management issues. When the computational modeling effort was initiated, this 
experiment was known as the Cryogenic Fluid Management Facility (ref. 1), 

(CFMF), but has since been renamed the Cryogenic On-Orbit Liquid Depot, Stor- 
age, Acquisition and Transfer Satellite (ref. 13), (COLD-SAT). The original 
design used a 0. 25-scale model of a propellant tank proposed by Boeing for a 
space based Orbit Transfer Vehicle, (OTV). This tank was selected as a proto- 
type for studying the effect of acceleration level on reorientation perform- 
ance. The design of the experiment relied on two Shuttle Reaction Control 
System (RCS) thrusters to provide an acceleration environment of 7.85 cm/s^, 
(8xl0"3g) . 

For the code validation phase, the investigator reviewed graphical dis- 
plays of the flow field evolution and judged settling to be achieved when a 
sufficient quantity of liquid had collected into a pool at the "bottom" of the 
tank. To eliminate the subjective nature of this evaluation, a measurable 
parameter was defined. Since the depth of the pool at the tank centerline can 
be tracked as the flow field evolves, the liquid is considered settled when 
this depth exceeds 20 percent of the total tank length. The settling time is 
defined as the elapsed time between initiation of thrust and satisfaction of 
this criterion. Although a few apparently anomakus settling times were pre- 
dicted at the lower thrust levels, the settling ciiterion worked well for the 
major! ty of cases . 

Figure 6 shows the shape and dimensions of the 0.25 tank as well as the 
computational mesh used to analyze the propellant motion within it. The analy- 
ses were performed for a tanx. 50 percent full of liquid hydrogen with ai' of 
the propellant initially collected at the top of the tank. Figure 6 shows the 
corresponding shape of the initial free surface. Figure 7 displays a sequence 
of flow fields which occur during the reorientation process with a SxlO'^g set- 
tling force. Tne liquid leading edge moves rapidly towards the bottom of the 
tank. A large geyser forms w’thin 8 sec of thrust initiation. The geysen is 
so severe that the liquid rebounds from the top of the tank before collecting 
into a pool after approximately !6 sec has elapsed. The liquid motion has 
been solviolent significant vapor pockets have been encapsulated in the pool. 

If this process orcirred in a real propulsion application, it might be nec- 
essary to extend the thruster firing period to insure that the vapor poci'ets 
are expelled from the pool before the main rocket is started. 

I. Sumner's rroort indicates that optimal ' eor i entation for the 0.25 tank 
should occur' for an acceleration level of appi o* Imakel y 3.6<10'^ cm/s*-. 
(3.7x10~^g). This level is ptedirted to be opkmal 'n the sense that it mini, 
mizes expenditure of propellant. Figure 8 presents a seqij?nce of flow fields 



m 





corresponding to this acceleration level. The liquid moves smoothly toward i| 

the bottom of the tank collecting into a sizeable pool within 1.5 min. A mod- 
erate geyser was formed but fewer vapor pockets were trapped within the pool. ' i 

After approximately 3 min. almost all of the liquid has collected in the bot- 
tom of the tank. I. Sumner (ref. 7) proposed that the reorientation process 
be judged complete when either the geyser settles back into the pool or the 
liquid film has cleared the tank wall. The propellant expenditure to accom- 
plish reorientation is roughly proportional to the vehicle delta-v Incurred 
during reorentat ion . The value to delta-v is easily computed by multiplying » 

the specified acceleration by the elapsed time required to accomplish reorien- 
tation. Using the RCS thrusters results in a vehicle delta-v of 125 cm/s ‘ 

whereas the optimal acceleration level corresponds to a delta-v of 6.5 cm/s, a *' 

fuel savings factor of almost 20:1. 

As the design of COLO-SAT evolved, the tank scale and shape were changed 
to more accurately emulate current OTV design concepts and to minimize tank , * 

thermal mass. The resulting tank geometry is presented in figure 9. This is 
a 0.215-scale model of a tank known as the Boeing Short SB OTV. Analyses 
therefore shifted to this new tank using the computational mesh shown in fig- 
ure 9. Acceleration environments between 2x10"^ and 1 cm/s^ were studied for ' 

a tank 50 percent full of liquid hydrogen. Figure 10 presents a sequence of f 

flow field depicting the reorientation process for an imposed acceleration of 
3.92x10"^ cm/s^, (4.00x'0^g). Settling time and vehicle delta-v were 
focused upon as the key parameters I'epresent i ng settling performance. 

f 

Settling time is of obvious interest for the scheduling of orbital maneu- V 

vers. Vehicle delta-v is used as measure of efficiency since it is directly 
correlated to the propellant expenditure. Figure 11 displays the relationship 
predicted between settling time and acceleration level. The anomolous set- 
tling times encountered at some of the lower accelerations have not yet been 
fully investigated, but are believed to be a resonance between the accelera- ' 

tion level and the geyser rise velocity. A somewhat arbitrary settling crite- 
ria was usea in this study and may also contribute to anomolous settling 
times. Figure 12 displays the relationship between vehicle delta-v and accel- 
eration level . | 

f 

FULL-SCALE OTV TANK | 

Upon completion of the small-scale tank cialyses, attention was focused on | 

a full-scale Boeing Short SB OTV. Modeling of the reor i entation processes in j 

this tank covered a range of acceleration environments from 1.57x10"'^ cm/s^ j 

(1.60xl0“'g) to 7.85x10"' cm/s^ (8.00xl0~‘* g) and included tanks 25, 50, and * 

75 percent filled with liquid hydrogen. The same mesh was used for these 
analyses as was used for the scale model of tne same shape. Figure 13 show: a 
sequence of velocity fields predicted for a 75 petcent full tank subject to an 
imposed acceleration environment of 1.96xl0"2 cm/s^ (2.00xl0"5 g). Figures '4 
to 16 display the relationship between settling time and acceleration level , 

for the three fill levels. Again, anomolies are encountered at the lower 
accelerations. Figures 17 to 19 display the relationship between vehicle 
delta-v and acceleration levei. The trends are not surprising, but ECLIPSE 
now provides a far more accurate tool for trading settling time versus propel- 
lant expenditure than was previously available. 
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DIMENSIONLESS SCALING OF REORIENTATION 


Since F.CLIPSE provides a tool capable of modeling the reor ientation proc- 
ess in both scale model tanks and full-scale tanks, it became possible to 
search for a dimensionless parameter capable of scaling experimental results 
from small-scale tanks to full-scale spacecraft tanks. In particular, the 
results for the scale model Boeing SB OTV tank and for the full-scale tank 
were used as the basis for this investigation. Following the investigators, 
the first attempt at correlation used the Bond number based on tank diameter 
as the scaling parameter. The data did not collapse using this scaling 
parameter . 

Various combinations of what were thought to be the relevant dimensionless 
variables (Bo, We, Fr, Re) were tried. The most successful attempt at corre- 
lating delta-V was a nondimensional grouping this paper defines as settling 
number . 


Se = Ra ) = Settling Number 


The Settling Number can be written as a function of more recognized dimen- 
sionless parameters: 

Se = (.^oWe > ’ 

'Re' ' ° 

When viewed in this lignt, it is seen as representing the ratio of viscous and 
gravitational forces to surface tension forces. 

The proof of a proposed correlating paramete’' is in the viewing the 
■'“suits. Figure 20 shows the relationshio between vehicle delta-/ and 5e for 
e Boeing 'jhort S3 OTV with a tank fil'Miq of 5'J ce'cert. A single straight 
line passes through all 24 data point,. t jnate i y . the analyses performed 

for the scale model tanx included only a few cases with 25 and 75 percent fill- 
ings. Although these analyses also correlate into a straight l^ne, they are 
too few to claim as support f^or the corielating parameter. They are however 
distinctly different lines from each other and fi'om the 50 percent case. 
Therefore, it appears tnat Se is a suitaple cor-elating paramete'- for relating 
reorientation performance from small-scale to full-scale geometrically similar 
tanks, containing the same fluid and with the same volume fraction of liquid. 


SUMMARr' OF RESULTS 


The ECLIPSE Code has been used to model the process of impulsive reorienta- 
tion. The accu^'acy of soniputationa 1 predictions was evaluated Pv compa''ison 
to e «per i men ta 1 1 y obtained data f()r i eor i enta t i on in smal'-sca’e tanks with 
shapes typical of spaceoaft PTopeilant tanks. The model correctlv predicted 
the extent of geyser fnmatior' and the elapsed time reqi,iired to accomplisn set- 
tling. Based on the c smpai i n rrs , ECLIPSE was C.idged ts> be a suitable tool for 
studying impulsive reor ient.at ion in ciysger'ic p' 0 [)e’-ant tanks. 

Reorientation of Hiuid hydiognn in flCght e<pe' 'ment tan^s was analyzed. 
These tanks have elliptical e.a d s ^ on n p i' « g g vy ( r v i 'ndf leal bar^'o! set^'on 
and are representat i ve jf ve-'ic'e pi pe''an- t-r-- k ;._ipgk was ap'p to mode' 
the reoi'ientation wi'-hoit di^fir.l*/ an,) p-o.iiej sif'i'ficant 'nsight into tne 
process. For one tank, i*" was dem<;r'"t'ate'J *k'ar if •■Mp e<’^ting des'g" was 


repldC8d with one producing an optimal acceleration environment propellant 
expenditure could be reduced by a factor of almost 20. For the other tank a 
range of acceleration environments was investigated and a summary of the 
results is reported. 

Reorientation in a full-scale OTV tank was modeled for three different 
tank fillings across a range of acceleration environments. A summary of the 
results for these cases is presented. 

A dimensionless parameter called the Settling number, Se . is proposed fo^ 
correlating the reorientation process between geometrically similar tanks with 
the same liquid volume fraction. To test the proposed parameter, computational 
predictions of vehicle delta-v acquired during settling for a full-scale space- 
craft tank and for a 0.2!5-scaIe model were plotted against Se. All data 
points fall into a single straight line, supporting the validity of Se as the 
appropriate correlating parameter. 
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TABLE I. - I.E. SUMNER'S TEST CONDITIONS MODELED USING ECLIPSE 
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